A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 
08(09.03.1,09.13.2,09.19.1,10.03.1,13.19.3) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Molecular gas in the Galactic center region 



oo 



III. Probing shocks in molecular cores* 

S. Huttemeister 1 ' 2 , G. Dahmen 3 ' 4 , R. Mauersberger 5 , C. Henkel 4 , T.L. Wilson 4 6 , J. Martin-Pintado 7 

1 Radioastronomisches Institut der Universitat Bonn, Auf dem Hiigel 71, D - 53121 Bonn, Germany 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, U.S.A. 

3 Physics Department, Queen Mary & Westfield College, University of London, Mile End Road, London El 4NS, England 

4 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D - 53121 Bonn, Germany 

5 Steward Observatory, The University of Arizona, Tucson, AZ 85721, U.S.A. 

6 Submillimeter Telescope Observatory, The University of Arizona, Tucson, AZ 85721, U.S.A. 

7 Centro Astronomico de Yebes, Apartado 148, 19080 Guadalajara, Spain 



in 



Received 1 Decemberl997 / Accepted 24 February 1998 



> 

o 

m 

o 

oo 

On 

S3 

Oh' 

6 

CO 



X 



Abstract. Multiline observations of C 18 and SiO iso- 
topomers toward 33 molecular peaks in the Galactic center 
region, taken at the SEST, JCMT and HHT telescopes, 
are presented. The C 18 presumably traces the total H2 
column density, while the SiO traces gas affected by shocks 
and high temperature chemistry. The J = 2 — > 1 line 
of SiO is seen only in few regions of the Galactic disk. 
This line is easily detected in all Galactic center sources 
observed. A comparison of the strength of the rare iso- 
topomers 29 SiO and 30 SiO to the strength of the main 
isotopomer 28 SiO implies that the J = 2 — > 1 transition of 
28 SiO is optically thick. The 29 Si/ 30 Si isotope ratio of 1.6 
in the Galactic center clouds is consistent with the terres- 
trial value. Large Velocity Gradient models show that the 
dense component (t*h 2 > 10 4 cm~ 3 ) in typical molecular 
cores in the Galactic center is cool (Tkin ~ 25 K), contrary 
to what is usually found in Giant Molecular Clouds in the 
disk, where the densest cores are the hottest. High kinetic 
temperatures, > 100 K, known to exist from NH 3 studies, 
are only present at lower gas densities of a few 10 3 cm~ 3 , 
where SiO is highly subthermally excited. Assuming that 
C 18 traces all of the molecular gas, it is found that in 
all cases but one, SiO emission is compatible with arising 
in gas at higher density that is (presently) relatively cool. 
The relative abundance of SiO is typically 10 -9 , but differs 
significantly between individual sources. It shows a depen- 
dence on the position of the source within the Galactic 
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center region. High abundances are found in those regions 
for which bar potential models predict a high likelihood 
for cloud-cloud collisions. These results can be used to re- 
late the amount of gas that has encountered shocks within 
the last ~ 10 6 years to the large scale kinematics in the 
inner ~ 500 pc of the Galaxy. 



Key words: ISM: clouds, molecules, structure - Galaxy: 
center - Radio Lines: ISM 

1. Introduction 

The molecular environment in the inner ~ 8° of the 
Galaxy differs drastically from that in the Galactic disk. 
Large scale surveys of CO isotopomers (e.g. Bally et al. 
1987, Heiligman 1987, Jackson et al. 1996, Dahmen et al. 
1997a (Paper I), Bitran et al. 1997) show that the cloud 
and the intercloud medium in the gaseous bulge of the 
Milky Way are molecular. The gas is characterized by 
large linewidths, indicating a high degree of turbulence. 
As shown by Dahmen et al. 1998 (Paper II), there is no 
simple relationship between CO line intensities and H2 
column densities and thus no global iV(H2)//(CO) con- 
version factor. 

Except for a few extraordinary regions such as Sgr A 
and SgrB2, little evidence for ongoing massive star forma- 
tion in Galactic center Giant Molecular Clouds (GMCs) is 
found, as is demonstrated by a general lack of strong FIR 
or radio continuum point sources associated with these 
clouds (Odenwald & Fazio 1984, Giisten 1989). Ambient 
dust temperatures are fairly low at Td ~ 25 — 30 K (e.g. 
Cox & Laurijs 1989). Based on ISO data, this is also a 
typical Td in the sources studied in this paper. 

In the Galactic disk, quasithermal SiO emission is 
tightly correlated with high temperature regions (e.g. 
Ziurys et al. 1989). A close association with outflows 



2 



S. Huttemeister et al.: Molecular Gas in the GC Region III. Probing shocked molecular gas 



strongly suggests that grain disruption by shocks is the 
major mechanism for releasing SiO into the gas phase 
(Martin-Pintado et al. 1992), although high temperature 
gas phase chemistry (Langer & Glassgold 1990) may play 
a minor role. In the Galactic center region, SiO is much 
more widespread (see e.g. the survey in the J = 1 — > 
transition by Martin-Pintado et al. 1997), which was in- 
terpreted as evidence for large scale or ubiquitous (fast) 
shocks. 

In this paper, C 18 and SiO data are presented for 
molecular cloud cores of the Bally et al. (1987) CS survey. 
Our measurements allow to (1) trace the H2 column den- 
sity, (2) determine or constrain density and temperature 
structure, (3) estimate SiO abundances and (4) obtain in- 
formation about silicon isotope ratios. Since our sources 
are selected on the basis of their intensity in CS, a gen- 
eral high density tracer, the sample is not a priori biased 
toward strong SiO emission. 

NH3 data (Huttemeister et al. 1993b) show that at 
least two phases of different kinetic temperatures are 
present within all cloud cores without massive star for- 
mation: A cool component with T^ n of 20-30K, close to 
the temperature of the dust in the Galactic center and a 
warm component with T^ n > 120 K. Tki n ~ 75 K is con- 
sidered by many authors, (e.g. the reviews by Morris & 
Serabyn 1996, Mezger et al. 1996) as 'typical'. This is, 
however, just the average over the hot and cool component 
and has no meaning as a distinct physical component. 

A major aim of this study is to determine the physical 
parameters, distribution and origin of these phases and to 
decide which of them, cool or warm, is associated with the 
bulk of the SiO emission. 



2. Observations 

The observations of the J = 1 — > and 2 — > 1 transitions 
of C 18 and the J = 2 -> 1 transition of 28 SiO, 29 SiO and 
30 SiO were carried out at the Swedish-ESO Submillimetre 
Telescope (SEST) at La Silla (Chile). All C 18 0(3 -> 2) 
data were taken at the James Clark Maxwell Telescope 
(JCMT) on Mauna Kea, (Hawaii), where we also obtained 
some 28 SiO(5 — > 4) and C 18 0(2 — > 1) measurements for 
calibration purposes. Finally, the 5^4 transition of 28 SiO 
was measured at the Heinrich Hertz Telescope (HHT) on 
Mt. Graham (Arizona). Table summarizes the observa- 
tions, including line frequencies, telescope parameters, re- 
ceivers and spectrometers used and typical system temper- 
atures. A telescope beamsize of 46" corresponds to 1.9 pc 
at a distance of the Galactic center of 8.5 kpc. 

Typically, the telescope pointing was checked every 
three hours. At SEST, two strong SiO masers close to the 
Galactic center, AHSco and VXSgr, were measured. At 
the JCMT, the pointing was done using five-point contin- 
uum maps on NGC 63341. At the HHT, continuum cross 
scans on Sgr B2 were used. At SEST and JCMT, the point- 



ing accuracy was always better than 5". At the HHT, the 
accuracy was always better than 10" and usually 5". 

At all telescopes, calibration was done with the stan- 
dard chopper wheel method, giving temperatures on the 
Tjj[ scale. These were converted to Tmb, the main beam 
brightness temperature (see, e.g., Downes 1989 and Rohlfs 
& Wilson 1996 for definitions), using the main beam effi- 
ciencies, ?7mb, listed in Table Q. 

To account for structure in our extended sources, we 
mapped the beam used to measure the C 18 0(1 — > 0) line 
in the (2 — > 1) and (3 — > 2) transitions. Due to time re- 
strictions, this was not possible for the SiO (5 — > 4) line 
observed at the HHT. From maps obtained at the JCMT 
for two sources in SiO (5 — > 4), the intensity in a 33" beam 
is higher by 8% and 15% than in 46" and 57" beams, re- 
spectively. This is caused by the structure of the individ- 
ual sources, and the numbers are, therefore, uncertain. We 
will, however, use these as correction factors when com- 
paring HHT SiO (5 — > 4) data to measurements obtained 
with larger beamsizes. 

To adjust the calibration scales of the three telescopes, 
we have observed a number of sources in the 1.3 mm 
band which is accessible to all instruments (see Table [j]). 
For C 18 0(2 — > 1), the integrated intensities at the 15m- 
telescopes SEST and JCMT agree to better than 10%, i.e. 
to within the uncertainty of the calibration of the indi- 
vidual measurement. The SiO(5 — > 4) JCMT maps, con- 
volved to the larger size of the HHT beam, also give in- 
tensities that are consistent to within 10% with the HHT 
data. Thus, the TMB-scales are in excellent agreement, 
and we are confident that our complete data set has been 
placed on a compatible intensity scale. 

We carried out all observations employing position 
switching. This was necessary given the large angular ex- 
tent of the sources, especially in C ls O. Offset positions 
were typically chosen 15' (SiO) or 30' (C ls O) away from 
the source, perpendicular to the Galactic plane. This ob- 
serving mode causes baselines that are not always per- 
fectly flat. Thus, polynomial baselines of order < 3 were 
subtracted from the data. Baselines were determined over 
a region of ±150 — 200 kms -1 from the line center. For the 
JCMT data, which were observed in exceptionally good 
weather (T225 < 0.05 throughout), linear baselines were 
usually sufficient. 

For the final analysis, all spectra were smoothed to a 
velocity resolution of ~ 2kms~ 1 . 

3. Results 

The 33 sources we study in detail are selected as dis- 
tinct peaks in the CS survey of Bally et al. (1987). They 
are, therefore, cores that should have large amounts of 
dense gas, which are embedded in a smooth, lower inten- 
sity molecular intercloud medium. Since we wish to derive 
the properties of typical clouds in the Galactic center envi- 
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Fig. 1. Sample spectra of eight sources with very complete observations. In the top panels, the C 18 lines are 
displayed: The (1 — > 0) transition is drawn as a dotted line, the (2 — ► 1) transition is drawn as a thin solid line and the 
(3 — ► 2) transition is displayed as a thick solid line. The bottom panels show the SiO spectra. 28 SiO(2 — > 1) is denoted 
as a thin solid line, 28 SiO(5 — > 4) is given by a thick solid line; the 29 SiO(2 — ► 1) spectrum is drawn as a thick dotted 
line and the 30 SiO(2 — > 1) is shown as a thin dotted line. The intensities of both rare isotopomers were multiplied by 
5 to allow clearer comparisons. Black squares overlayed on the map of the overall distribution of the J = 1 — > line 
of C 18 (taken from Paper I) mark all 33 sources we observed except Clump 1 at I ~ —5°. 
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Table 1. Summary of the observations 
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a) : SSB: Single sideband receiver, with a sideband rejection of ~ 20 dB ; DSB: Double sideband receiver 

b) : on a Tmb scale 

c) : AOS: Acousto-Optical Spectrometer; DAS: Dutch Autocorrelation Spectrometer 

d) : Only three sources were observed for cross-telescope calibration with SEST 

e) : Only two sources were observed for cross-telescope calibration with the HHT 



ronment, the non- typical molecular peaks associated with 
Sgr A and SgrB2 were not included in the sample. 

The (1 -> 0) and (2 -> 1) transitions of C 18 and 
the (2 — * 1) transition of 28 SiO were observed toward all 
sources of our sample. Both the C 18 and (more notably) 
the SiO lines were easily seen in all sources. We observed 
the 10 strongest C 18 sources in the (3 — > 2) submil- 
limeter transition of C 18 0, again detecting all. The rare 
isotopomers 29 SiO and 30 SiO were observed and detected 
toward 12 clouds showing strong emission in 28 SiO. Of 
8 sources measured in the 28 SiO(5 — > 4) line, 7 showed 
emission. In Fig.0 and Table|, we give an overview of 
our results. Sample spectra for eight sources are displayed 
in Fig.]!]. We mark the position of all sources observed 
on a survey map of the (1 — ► 0) transition of C 18 (Pa- 
pers I and II). Galactic coordinates are used for source 
names. For the corresponding equatorial coordinates, re- 
fer to Huttemeister et al. (1993b). Most of the clouds we 
have studied are located within the molecular bulge re- 
gion, while three lie within the 'Clump 2' complex (Bania 
1977, Stark & Bania 1986) at I ~ 3°. 'Clump 1' is located 
~ 5° south of the Galactic center. 

Integrated line intensities for the C 18 0(2 — > 1) and 
28 SiO(2 — > 1) transitions and line intensity ratios derived 
from integrated intensities are presented in Table |[ The 
errors are determined from the (usually small) formal er- 
ror of a Gaussian fit to the lines or, in case of non-Gaussian 
lincshapes, from the rms noise in the spectrum. In addi- 
tion, a calibration uncertainty of 10% for data taken with 
the same telescope or 20% for data obtained with different 
instruments was assumed. For Gaussian lines, the formal 
error of the fit agrees closely with the error obtained from 
the rms noise in the spectrum. 



For a given species, the center velocities and line- 
shapes of the different transitions always agree, to within 
the noise. Between SiO and C 18 0, however, there can 
be significant differences. The sources M+0.50+0.00 and 
M+0.24+0.02 are good examples: In both cases, the C 18 
lines show two distinct peaks, while the SiO transitions are 
single-peaked. In M+0.50+0.00, the central velocity of the 
SiO is close to the weaker C 18 line; in M+0.24+0.02 SiO 
and the stronger C 18 peak agree. It is also noteworthy 
that narrow C 18 lines close to a wlsr of Okms" 1 , likely 
of local origin, never have a counterpart in SiO. This is il- 
lustrated by the source M-0.51-0.16 in Fig.|l| and demon- 
strates the unusual nature of the Galactic center sources 
as compared to Galactic disk clouds. 

Since the lines of 29 SiO and 30 SiO are very likely op- 
tically thin, we can directly check whether the ratio of 
29 Si/ 30 Si in the Galactic center region agrees with the 
terrestrial value. We find a line intensity ratio of 1.6 ± 0.2 
in our sample, in excellent agreement with the terrestrial 
isotope ratio of 1.5. This confirms that this ratio does not 
depend on the galactocentric distance (Wilson & Rood 
1994, Penzias 1981), which is an expected result if both 
isotopes are synthezised in the s-process in stars of the 
same type. 

For 28 SiO/ 29 SiO, we take the ratio to be the terres- 
trial value of 20, as suggested by Penzias (1981). Since 
the 28 SiO(2 — ► 1) transition is optically thick, we cannot 
check this assumption, but it is supported by the largest 
line ratio pfg9 2 g (toward the exceptional source M+1.31- 
0.13, see Section 4.3), which is indeed close to 20. 



4. The physical conditions of the gas 
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Table 2. Integrated intensities of the C 18 0(2 — > 1) and SiO(2 — > 1) transitions and line intensity ratios of C 18 and 
SiO. 
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Abbreviations used throughout the paper: 7^°= 7 ( c 2 _^ ^ /2 °= /g_ 2 ^// ( C 2 _^, ^ 5/ ^ i0 = / (B 5°/J(2-i° 

„SiO _ r 28 SiO /r 29 SiO „SiO _ r 29 SiO / r 30 SiO 
' 28/29 — (2 — 1) / (2 — 1) ' ' 29/30 — J (2-l) / J (2-l) 

a) : The result is uncertain, due to either baseline uncertainties or low signal-to-noise ratios 

b) : Non-gaussian lincshape or several blended velocity components; ratios have been calculated by adding channels with 
emission 

c) : very narrow line, possibly from local ISM 

d) : 3cr 
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4--1- Determinination of densities and temperatures 

A Large Velocity Gradient (hereafter LVG) model was 
used to estimate physical conditions of the gas, namely 
density (n(H2)), kinetic temperature (Tki n ) and - in the 
case of SiO - molecular abundance (AT(SiO)) (see Scoville 
& Solomon 1974 and De Jong et al. 1975 for the general 
properties). We used a modified version of the code of 
Henkel et al. (1980). The basic LVG assumption is that 
a of systematic velocity gradient large compared to ran- 
dom motions and monotonic. This allows us to treat the 
molecular excitation as a local problem. This is certainly 
an idealisation for Galactic center clouds. However, an ap- 
plication of an LVG code requires no detailed knowledge 
of the velocity field. In the case of CO lines, H2 densities 
determined under the assumption of LVG and microtur- 
bulence (the opposite extreme) do not differ by more than 
a factor of three (White 1977). 

Our analysis is also based on the assumption that, for a 
given molecule, all transitions observed arise in the same 
volume. Since the excitation depends on critical density 
and optical depth, and thus is not identical for different 
transitions or different isotopomers, this assumption may 
not be strictly correct, especially if the medium is very 
clumpy (see, e.g., the discussion in Oka et al. 1997 for 
12 CO). For C 18 0, this is not a critical problem since the 
LVG calculations show that all transitions measured by us 
are optically thin, or, in the low temperature/high density 
scenario described below, reach, at most, r w 1. There- 
fore, C 18 line intensity ratios are not strongly affected 
by the possibly different locations of cloud 'photospheres' 
from which the bulk of the photons are emitted. All tran- 
sitions of C 18 should trace all molecular gas at densities 
> 10 3 cm~ 3 , as long as C 18 is not selectively dissoci- 
ated by UV radiation. Away from the star forming regions 
SgrA and SgrB2, the UV radiation field in the Galactic 
center region is not likely to be strong (e.g. Nagakawa et 
al. 1995). In addition, because we are considering density 
peaks of the molecular gas (and dust), the destruction of 
C 18 should not be a problem; from the H2 column den- 
sities derived below, dust extinction is likely to provide 
sufficient shielding. 

In our model calculations, the population n(J), exci- 
tation temperature T ex , brightness temperature Tb, and 
optical depth r, are determined for each rotational level up 
to a maximum J = 20 for C 18 and J — 13 for SiO. Col- 
lision rates are taken from Green & Chapman (1978). We 
have chosen a constant velocity gradient of 5kms _1 pc _1 , 
a value estimated from the velocity extent and diameter of 
Galactic center clouds as mapped by Bally et al. (1987) or 
Lindquist et al. (1995). The background continuum tem- 
perature was assumed to be 2.7 K. 
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Fig. 2. Line intensity ratios from our LVG model calcula- 
tions for C 18 (upper panel) and 28 SiO and 29 SiO (lower 
panel). The contours are at: 

llf/l (solid): 0.3, 0.45, 0.6, 1.0, 1.5, 2.0, 3.0, 3.5, 3.75 

nf^° (dashed): 0.1, 0.2, 0.35, 0.5, 0.75, 1.0, 1.5, 2.0 

ftg/f (solid): 0.04, 0.05, 0.1, 0.2, 0.5, 0.8, 0.9, 1.0 

r 28/2g( dashed ) : L2 > L5 > 2 -°> 5 -°> 10 -°> 15 -°> 19 -°> 19 - 75 > 
20.0 

4.1.1. C 18 

We assume that the abundance of C 18 relative to H2 
is 4 • 10 -7 . This ratio is based on the standard Galactic 
center 16 0/ 18 isotopic ratio of 250 (Wilson & Matteucci 
1992, Wilson & Rood 1994) and a CO/H 2 ratio of 10™ 4 
(e.g. Frerking et al. 1982, Duvert et al. 1986 (dark clouds), 
Blake et al. 1987 (OMC-1)). According to the calculations 
of Farquhar et al. (1994), this ratio is stable against a 
possibly enhanced cosmic ray flux close to the Galactic 
center. For C 18 0, we then calculate a grid by varying the 
H2 density, n(H2), from 10 2 cm~ 3 to 10 6 cm" 3 and the 
kinetic temperature, Tki n , from 5 to 200 K. 

We found the range of n (H2) and T^ n by comparing 
the observed line intensity ratios (Table|^) to the ratios 
predicted by the model (Fig.||). It is important to use line 
intensity ratios, not just intensities, since intensities are 
affected by an a priori unknown beam filling factor, /b. 
For the intensity ratios, we took /b to be identical for all 
transitions of a given molecule. 
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The integrated intensity ratio of the C 18 0(2 — > 1) to 
(1 — ► 0) transition, 7^2/1°' f rom au sources is 1.3 ± 0.5, 
with extreme values ranging from 0.5 to ~ 2. The ra- 
tio of the (3 -> 2) to the (2 -> 1) intensity, llfl , is 
0.5 ± 0.2. Averaged over all sources, the LVG model fits 
yield T cx = 11 (±4) K for the C 18 0(2 -> 1) transition, that 
is, C 18 is subthermally excited. The beam averaged peak 
optical depth T21 is always < 0.2. The measured peak 
temperatures Tmb indicate a beam filling factor 
of - 1 for n(H 2 )= (2 - 4)10 3 c m - 3 . fg 18 ° is lower if 
n (H 2 ) is higher. From the structure seen in our maps of 
the C 18 0(1— >0) beam in higher transitions, we estimate 
/b °, averaged over all gas components, to be ~ 0.5. 
These maps, however, have a low signal-to-noise ratio. 

From the model calculations (Fig.|2|, upper panel), 
maximum values up to 4 in 7^2/1° ano - U P ^° ^ m ^3/2° 
are expected for hot (Tk; n > 100 K) and dense (n(H2)~ 
10 4 — 10 5 cm -3 ) gas. Such high ratios are reached in none 
of our sources. Possible combinations of n (H 2 ) and Tk; n 
are either a high Tki n of > 100 K and n (H2) of (1 - 
4) 10 3 cm" 3 or a lower T kin (< 20 K - 30 K) at a higher 
^(H^) > 10 4 cm~ 3 . We have carried out this analysis in- 
dividually for all clouds. For all sources our C 18 data 
unambiguously require gas with densities of > 10 4 cm~ 3 
in these Galactic center clouds to be cool. 

From ammonia measurements (Hiittemeister et al. 
1993b) it is known that both hot and cool gas is present 
within the area covered by our beam. Thus, it is likely that 
both scenarios possible from the C 18 data are realized. 
A continuum of temperatures and densities is compatible 
with our data. Then Tki n rises steadily as n (H2) decreases, 
when the line of sight samples different parts of the cloud. 
Density peaks (indicated by strong CS lines) in Galactic 
center region GMCs not presently undergoing massive star 
formation can be regarded as 'cool dense cores' 'in contrast 
to the hot dense cores generally encountered within GMCs 
in the disk. 

4.1.2. SiO 

It has been shown that the fractional abundance of SiO, 
A (SiO), can change by more than six orders of magnitude 
between quiescent cold material and hot, shocked dense 
gas (Ziurys et al. 1989, Martin-Pintado et al. 1992). Since 
SiO is readily detected in all our sources, A(SiO) cannot 
be low. However, we cannot assume any specific value. 
Therefore, we have used A(SiO) as a free parameter in 
the models, varying it from 10~ 12 to 10~ 6 . Grids have 
been calculated as a function of n.(H2) and A(SiO) for 
T kin = 25 K, 50 K, 75 K, 100 K, 150 K. 

The average line intensity ratio between the SiO (5 — * 
4) and (2 -> 1) transition is TlZf = 0.16±0.05, for those 
sources where we detected SiO(5 — ► 4). M+1.31-Q.13 is ex- 
ceptional and will be discussed further in Section 4.3. Since 
the signal-to-noise (S/N) ratio in the (5 — » 4) line is not 



always high and the (2 — » 1) SEST beam was not mapped 
in this transition, it is not clear whether source to source 
differences (apart from M+l. 31-013) are significant. For 
the line ratios between the (2 — > 1) transitions of 28 SiO 
and 29 SiO (excluding M+l. 31-0. 13), we find a mean ratio 
of rfgj 2g = 12.2 ± 2.2 with extreme values ranging from 
7.4 to 15.9. Since these data, even in the 29 SiO transition, 
have very good S/N ratios and were measured with the 
same telescope and beamsize, there is no doubt that the 
source-to-source differences are significant. 

Fig. [2] (lower panel) shows the line ratios T^ 5 ^° and 

r 28/29 f° r ^kin= 75 K from our LVG calculations. Lowering 
Tkin to 25 K requires a slight increase of the H2 densities 
for a given line ratio, by a factor of ~ 1.25. Raising Tki n to 
150 K corresponds to slightly lower densities. In general, 
the effects of changing Tkin on the line ratios are negligible: 
SiO line ratios are almost insensitive to cloud temperature. 

To derive T ox and n (H 2 ), we have analysed all sources 
individually. Excluding M+l. 31-0. 13, we find T ex ranging 
from 5K to 10 K. n(H 2 ) changes with the assumed SiO 
abundance: A higher A(SiO) requires a lower H2 density 
(Fig.||). Curves of constant TZ 5 ^ l ° and r|g9 2g run paral- 
lel for a range of densities and SiO abundances. Typical 
ranges are A(SiO) « 10~ 6 for n (H 2 ) w (1 - 2.5) 10 3 cm" 3 
to A(SiO)w 10~ 9 for n(H 2 ) « 410 4 - 10 5 c m - 3 . Lower 
SiO abundances or higher densities are not possible, since 
then the curves diverge. 

Plotting the measured T M b's for the 28 SiO(2 -> 1) 
transition in the LVG plots (e.g. Fig.||, lower left panel), 
we find that a considerably higher Tb is needed to fit the 
observed line ratios. Thus we estimate that the beam fill- 
ing factor /f i0 (T MB = /b i0 r B ) is ~ 0.2 for the 'typical' 
sources. 

This analysis is biased toward more intense SiO lines, 
since only those were observed in transitions other than 
28 SiO(2 — > 1). The strength in SiO toward these positions 
may be either due to a large amount of molecular gas, or to 
an abundance of SiO that is above average. This question 
will be addressed in the next section. 

4-2. A joint view of C ls O and SiO 

4.2.1. Column densities and SiO abundances 

Further insight in the structure of the clouds is gained by 
combining what can be learned from C 18 and SiO. 

We calculate the total (beam averaged) column densi- 
ties for C 18 and 28 SiO, applying the same procedure to 
both molecules: First, we use the observed Tmb and the 
best fit T GX from the LVG model for the (2 — > 1) transi- 
tions (the lines with the best S/N ratio) to calculate the 
(beam averaged) optical depth f . Taking J f(v)dv = 
1.06 fo Avi/2, we then derive the column density in the 
J = 1 level, N\. The total beam averaged column density 
N is determined by dividing N± by the fraction of the pop- 
ulation residing in the J = 1 level, derived from the LVG 
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model. (See Rohlfs & Wilson (1996) for a collection of the 
(standard) formulae we used.) For non-gaussian lines, the 
second moment of the line has been used as a measure of 
the line width instead of Av 1 / 2 . This has been combined 
with a Tmb which reproduces the observed integrated in- 
tensity. This is correct as long as T cx does not change 
across the line profile. Since the line shapes in all transi- 
tions observed for a given molecule agree closely, this is a 
reasonable assumption. 

For the sources where we observed only the SiO(2 — > 1) 
line, we have used average values for T ex and the fraction 
of the population in the J = 1 level. These are 7 K and 
0.36, respectively. Total beam averaged column densities 
are given in Table |^. 

The usually optically thin, easily excited transitions 
of C 18 are known to be excellent tracers of H2 column 
density, 7Vh 2 , over a wide variety of H2 densities and ki- 
netic temperatures. In particular, the integrated intensity 
of the J = 2 — > 1 transition can be directly related to 
Nn 2 (Genzel 1992). Our LVG modelling confirms that for 
a large range of temperatures the fraction of the total pop- 
ulation in the C 18 J = 1 level remains almost constant 
between n (H2) = 10 2 cm~ 3 and 10 4 cm -3 . 

In most cases, we find that column densities agree with 
the formula given by Mauersberger et al. (1992) and Gen- 
zel (1992) (iV H2 = 1.12 10 21 /rCg O (2 -> l)dv for the 
Galactic center isotopic ratios) to within 10%. Differences 
up to a factor of 2 are found if the excitation temperature 
is exceptionally low. 

The knowledge of Nn 2 and N Sl ° allows us to derive 
a beam averaged SiO relative abundance, assuming that 
the C 18 and SiO emission arises from the same gas. This 
is given in Table|[ and typically ranges from 0.5 10 -9 to 
5 10 -9 . Since C 18 presumably traces allH2 with n (H2) > 
10 3 cm -3 , the true abundance of SiO cannot be lower than 
the value thus determined. It can be higher, if only part 
of the gas within the beam is abundant in SiO, which, 
considering the differences in line centers and line profiles 
between C 18 and SiO (e.g. M+0.24+0.02, see Fig.§, is 
likely. 

For those sources where the SiO(2 — » 1) transition is 
weak and which were only observed in this SiO transition, 
the H2 column density obtained from C 18 is not system- 
atically lower than for the more intense sources. On aver- 
age, in these clouds X(SiO) is lower by a factor of ~3. 
Source-to-source variations are significant, even though 
the properties in dense gas as traced by CS (Huttemeister 
1993) and total gas as traced by C 18 do not differ. This 
agrees with the variation found in the survey of Martin- 
Pintado et al. (1997). Note, however, that we find SiO 
with an abundance that is high when compared to typical 
disk GMCs in all sources. 



Table 3. Total beam averaged column densities of 
C 18 (iV cl8 °) and SiO (N sio ) and beam averaged abun- 
dances of SiO, X (SiO). iV H2 follows from N c °, assuming 
X(C 18 0) to be 4 • 10~ 7 . Errors are given in parentheses. 
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4.2.2. Density and temperature 

We now combine the LVG models for C 18 and SiO to 
derive a consistent solution for n(H2) and 7ki n . This is 
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done in terms of a two component model, where part of 
each C 18 emission line is ascribed to a cool component, 
with the remaining part arising from a hot component. 

From the NH3 studies of Hiittemeister et al. (1993b), 
we know that, for a typical cloud, roughly 25% of the 
neutral gas mass is contained in the hot (> 100 K) com- 
ponent, which, from our C 18 analysis, must have a low 
density. Thus, the remaining ~ 75% of the total gas mass 
is at a density of > 10 4 cm" 3 . The hot, thin and the cool, 
dense component are roughly in pressure equilibrium. Let 
us assume that the SiO emission arises in only the cool, 
dense component. The corresponding relative abundances 
of SiO, X(SiO), are 25% above the beam averaged values 
given in Table|[ X(SiO) lies within the range allowed by 
the SiO line ratios (Fig. ||) and is consistent with an H2 
density of typically (1 - 4) 10 4 cm" 3 . The beam filling fac- 
tor /b in this dense component must be the same for both 
the SiO and the part of the C 18 emission that arises in 
this component. 

We can estimate the main beam brightness temper- 
ature of the dense, C 18 emitting gas to be 0.75 of 
the total C 18 intensity. Using the beam filling fac- 
tor determined for the dense gas from SiO lines, the 
brightness temperature of the C 18 0(2 — > 1) line is 
T^ 18 °(75%)«0.75T MB //b- This simple relation can be 
used since the C 18 0(2 — > 1) line has an optical depth 
< 1. If the solution is self-consistent, the H2 density in- 
dicated by this value of Tb and the observed line ratios 
must agree with the density derived for SiO. To within 
a factor of 2, we find that this is indeed the case. Thus, 
it is possible that the SiO emission originates in gas hav- 
ing n (H 2 ) ~ (1 — 4)10 4 cm" 3 . This gas must be relatively 
cool, which is a surprising result in the light of the clear 
association of the SiO emission with hot shocked gas in 
the Galactic disk. The intrinsic optical depth tq of the 
SiO(2 — > 1) then ranges from 4-8. 

Can the SiO emission be associated with the hot, 
thin gas component? To check this, we used the same 
procedure as before: From the assumption that the SiO 
emission arises in ~ 25% of the gas, we get a value for 
X(SiO) and a corresponding H2 density. The beam fill- 
ing factor /b for SiO and, consequently, for the part of 
C 18 originating in the same gas does not change. Now, 
T B jl8 °(25%)« 0.25T MB /b- The additional requirement 
that this gas should be hot yields a H2 density that is, 
for many clouds, lower by an order of magnitude than 
what is derived for the same component from SiO. Thus, 
this is not a consistent solution. If the SiO arises in a still 
smaller portion of the total gas, the discrepancy becomes 
even worse. Note that we do not claim that the hot, thin 
component is devoid of SiO. We consider it likely that the 
SiO abundance in this component is similar to what we de- 
termine for the dense component. However, the SiO exci- 
tation in this component is extremely subthermal, leading 
(1) to lower line intensities than those we observe (specif- 



ically, no emission from the (5— >4) transition) and (2) to 
different line intensity ratios. 

Finally, we can rule out the presence of a hot, high 
density (ra(H2)> 10 4 cm" 3 ) component as the origin of 
the SiO emission: We know from the line ratios T^/f an d 
r 28/29 t na ^ -X"(SiO) must decrease with increasing density, 
independent of Tki n . Thus, a hot, high density component 
containing the bulk of the SiO would have to contain a 
significant fraction of the total H2 column density and 
mass. Such a component should be visible in C 18 0. Thus 
we conclude that it does not exist in these clouds. 

Therefore the only consistent scenario that accounts 
for both the SiO and C 18 results requires that most of 
the SiO emission arises from a cool (20 K - 30 K) gas com- 
ponent with n (H2) ~ 10 4 cm" 3 . 

So far, we have assumed that the C 18 abundance in 
all components considered is known and constant. If this is 
not the case, i.e. if C 18 is not a reliable H2 column den- 
sity tracer, the line of argumentation given above cannot 
be maintained. In this case, it becomes possible to claim 
that the SiO emission arises in a hot, dense component 
not seen in C 18 0. This is, however, only possible if C 18 
is selectively underabundant in the component where the 
SiO abundance is high. Since selective C 18 dissociation 
due to lack of shielding is efficient only in diffuse, low den- 
sity gas (if the UV field is strong enough), it is difficult 
to find a scenario causing C 18 to be underabundant in 
dense, warm gas. 

4-3. A detailed investigation of two sources 

M+0.83-0.18 In Fig. | (left panels), we illustrate the pro- 
cedure of deriving the physical cloud parameters described 
above for the source M+0.83-0.18. This is representa- 
tive of the type of cloud where SiO emission arises in a 
cool, moderately dense component. The SiO line ratios 
in this source (lower left panel in Fig.||) are well fit by 
T cx = 10K. Then, X(SiO) is 3.9 10" 9 if SiO arises from 
75% of the total gas. The observed T M b is 0.48 K, while 
a Tb of 7.5 K is required to fit the range of n(H 2 ) and 
X(SiO). This results in the lowest value for /b of all 
sources, 0.06, and the highest optical depth r , 8. Applying 
the same /b to the C 18 data (upper left panel), we have 
T b (C 18 O)=0.75T M b/0.06 = 5.1K. If we take the result- 
ing n (H2) from the SiO analysis, 4 10 4 cm -3 , and transfer 
this to C 18 0, we find that Tki n is only ~ 12 K. If we re- 
quire Tkin to be identical to the value we have measured for 
the cool component from NH3, 21 K (Hiittemeister et al. 
1993b), the density obtained for C ls O is 2 10 4 cm" 3 . Con- 
sidering all the assumptions, there is satisfactory agree- 
ment. Could the results be consistent with the 'high tem- 
perature scenario' described in the last section? We find 
that n (H2) required from an analysis of the SiO data is 
~ 3 10 4 cm" 3 , while the C 18 data at the appropriate Tb 
yield ~ 3 10 3 cm -3 , clearly a far less consistent solution. 
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Fig. 3. LVG Model calculations for two individual sources, demonstrating the range of physical conditions encountered. 
The observed line intensity ratios are plotted including their 1 a uncertainties. The Tb- curves are the observed Tmb 
and the best fit temperature derived from the line ratios. The coresponding beam filling factor /b is also given. For 
M+l.31-0.13, the observed and the best fit temperature are the same, thus /b = 1. t 21 is the peak optical depth of 
the 28 SiO line (not beam averaged). 



M+l.31-0.13 Now we analyze the one source which is 
very different from all the others in our sample: M+l.31- 
0.13 (right panels of Fig. [3]). This is the only source where 
the (5 — » 4) transition of 28 SiO was not detected, even 
though the (2 — > 1) line is strong. It also is the one source 
where rfg9 2 g approaches the terrestrial isotopic ratio of 
20, indicating that ro for 28 SiO(2 — » 1) cannot be large. 
M+l.31-0.13 has two velocity components. We analyze 
the (slightly) more intense wlsr ~ 85kms _1 component, 
but the «lsr ~ 35 km s -1 component shows the same char- 
acteristics. 

From the SiO line ratios, we derive an extremely low 
excitation temperature, 3.2 K. The large column density 
for SiO depends strongly on the exact value of T ex . There- 
fore, we can give only an estimate of 7V Sl ° in Table [|. The 
corresponding X(SiO) is ~ 10~ 8 , by far the highest in 
the sample. In addition, for this source we find Tmb = Tb, 
i.e. /b = I. As expected, this results in a lower r , ~ 2, 
in the 28 SiO(2 — > 1) line. For the H 2 density, we find 



n (Ha) ~ 2 10 3 cm 3 . This is a factor of ~ 10 lower than 
for all other sources. 

Since our calculations give r ~ 2 and /b = 1 for SiO, 
we conclude that in this source SiO traces all the gas, 
as does C 18 0. Hence, /b = 1 and Tb =Tmb is also valid 
for C 18 0. The resulting n (H 2 ) from C 18 is 210 3 cm- 3 , 
identical to the result obtained from SiO. As expected, 
Tkm is high, > 100 K. Thus, for this one source, we have 
found a fully self-consistent solution of an entirely different 
type: The SiO arises in thin, hot gas and has an abundance 
that is far above average, even for the clouds that are rich 
in SiO. 

Is there evidence that M+l.31-0.13 lacks the dense 
and cool gas component? Remarkably, our NH3 data show 
that this cloud is the only one in the entire sample in 
which the (J,K) = (4,4) inversion transition at 200 K 
above ground is more intense than the (1,1) transition at 
23 K above ground. Thus, in this cloud a cool, dense gas 
component containing the bulk of the gas is not present. 
Most, perhaps all of the gas is hot and thin. Relating this 
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to the high SiO abundance, it is likely that SiO does indeed 
form at high temperatures in the Galactic center region. 
In M+l.31-0.13, the SiO formation process is either still 
ongoing or has occured very recently. The SiO has not yet 
had time to recondense onto dust grains, and the cloud 
itself has not had time to form dense cool cores. 

5. Discussion 

5.1. Cloud properties as a tracer of large scale dynamics 

We will now explore how the relative SiO abundances may 
be related to the position of the cloud in the large scale 
Galactic center environment. 

In Fig.0 we plot the the beam averaged SiO abun- 
dance, X s , on a longitude-velocity diagram of the large 
scale distribution of C 18 0. While any estimates of line-of- 
sight locations in the Galactic center are rather uncertain, 
this plot gives a measure of location. 

There is a clear trend in the large-scale distribution of 
SiO abundances. For — 0?35 < I < 0?6, between the loca- 
tion of the Sgr C and Sgr B2 regions, we find a pronounced 
lack of very large A (SiO). Most clouds with very high 
abundances, including the exceptional source M+l.31- 
0.13, are located at I > 0?8. The 'Clump 2' Region is not 
included in this plot, since, at i » 3°, this is not really 
part of the continuous molecular bulge or, in the nota- 
tion of Morris & Serabyn (1996), the Central Molecular 
Zone (CMZ), although it shares many characteristics of 
the Galactic center gas. X Sl ° in this region ranges from 
1 10 -9 to 2 1CP 9 , typical 'Galactic center' values that are 
high compared to the disk but not as extreme as what is 
found at 0?8 < I < 1?5. 

The detection of an 'SiO hole' at -0?35 < I < 0?6 is 
confirmed by an extension of the survey of Martin-Pintado 
et al. (1997) (unpublished data). 

In recent years, the large scale dynamics of the gas 
in the CMZ, characterized by large non-circular motions 
and a distinct 'parallelogram' shape of ^-diagrams based 
on 12 CO and 13 CO (e.g. Bally et al. 1988), have been ex- 
plained in terms of a model involving a rotating (stellar) 
bar with corotation at 2.4 kpc and oriented at an angle of 
~ 20° with the line of sight to the Galactic center (see 
Morris & Serabyn 1996 and references therein). Within 
such a potential, there exist two categories of closed ellip- 
tical orbits, called x\ and a;2-type orbits. Inside a cusped 
orbit, the xi-orbits, which are elongated along the bar 
axis, become self-intersecting. Clouds on these orbits en- 
counter a shock and within a dynamical time plunge into 
£2-orbits that lie considerably deeper within the potential 
and mimic circular orbits (Mulder & Liem 1986, Athanas- 
soula 1988, 1992). This encounter breaks the flow along 
the cusped orbit into a spray that spreads out into the in- 
terior of the orbit (Binney et al. 1991, Athanassoula 1992, 
Jenkins & Binney 1994). 

According to this model, Sgr B2 and Sgr C can be inter- 
preted to be at the locations of the intersections of the x\- 
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Fig. 4. The beam averaged SiO abundance plotted on a 
longitude- velocity diagram of the large scale distribution 
of C 18 0, adapted from Paper I (spatial resolution 9'). 

and the X2-orbits. The dense, virialized clouds on the x 2 - 
orbits, inside of SgrB2 and SgrC, are more likely to form 
stars, while the less dense clouds on the x\ orbits are too 
disrupted by shocks for efficient star formation. Generally, 
this is supported by the fact that stars of young and in- 
termediate age are restricted to the area inside of Sgr B2 
and SgrC, suggesting sustained star formation in this re- 
gion (Serabyn & Morris 1996). In SgrB2, the ongoing star 
formation, the large amount of molecular gas, the exis- 
tence of a very hot molecular component (Hiittemeister 
et al. 1993a, 1995, Flower et al. 1995) and the high op- 
tical depth of the 12 CO emission (Paper II) also provide 
ample evidence for this picture. The Sgr C complex is a 
much more quiescent region, but it contains a large HII 
region and GMCs. 

Because the sprayed gas crashes into material that is 
still on xi-orbits, collision regions are expected along the 
acceleration part of the cusped orbits. Due to perspective, 
these main collision areas should be at higher positive lon- 
gitudes than SgrB2 and between Sgr A and SgrC. While 
the situation at negative I is confused (the gas on x 2 orbits 
and the collision regions are expected at similar Galac- 
tic longitudes), at positive I, especially the T?5-complex' 
(Bally et al. 1988) shows all characteristics expected from 
a collision region (see the discussion in Paper II). 

Inspecting Fig.|^ with this picture in mind, we find 
remarkable support from our data. The exceptional cloud 
M+l.31-0.13 and three other sources with high SiO abun- 
dances are located within the 1.5°-complex, exactly where 
strong shocks are expected. Generally, the high SiO abun- 
dances at I > 0?8 can be identified with a region in which 
the probability of gas encountering shocks is high. Thus, 
for the first time, we can tie the physical and chemical pa- 
rameters of Galactic center molecular peaks to the large 
scale dynamics of the region. 
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Since CS traces all dense gas, not just the part that 
has been subjected to shocks, we do not expect a similar 
'zone of avoidance'. Indeed, a plot similar to Fig [I| compar- 
ing the CS abundance in these clouds (data presented in 
Hiittemeister 1993), to the total H2 column density does 
not show such an effect. 

Of course, the correspondence is not absolute. Not ev- 
ery cloud at I > 0?8 must experience strong shocks. Note, 
however, that the clouds with very low SiO abundances in 
this range of I have negative velocities and might not be 
located in the collision region. Also, there are certainly 
strong shocks occuring in the Galactic center that are 
unrelated to large scale dynamics. Martin-Pintado et al. 
(1997) explain the 'SiO clouds' they find by shocks of a 
variety of origins: Interaction with SNRs close to Sgr A, 
interaction with non-thermal filaments in the radio arc 
and cloud-cloud collisions (as expected from large scale 
dynamics) or expanding bubbles in the vicinity of SgrB2. 
Their map extends from — 0?2 to +0?8. Therefore, they 
have missed the systematic signature of large scale effects 
present at higher positive longitudes. 

Another mechanism explaining the preferential oc- 
curence of shocks in the l?5-complex might be fossil su- 
perbubbles, remnants of a phase of Sgr B2 type star forma- 
tion activity. The expansion of superbubbles would offer a 
natural explanation of the large extent of the CO emission 
to positive Galactic latitudes perpendicular to the Glactic 
plane seen in this region. However, while the bar model as 
outlined above is two-dimensional and thus does not natu- 
rally produce vertical structure, sprayed gas colliding with 
gas on a; 1 -orbits might give rise to turbulence also push- 
ing gas out of the plane. It is noteworthy that a 12 CO 
map (see, e.g., Paper II and Bitran et al. 1997) shows the 
molecular gas extending toward negative b at the southern 
edge of the CMZ. Thus, it seems possible that the CMZ 
is warped or susceptible to instabilities close to its edges 
(Morris & Serabyn 1996 and references therein). 

The 'bar model' is not the only possibility to explain 
the large scale dynamics of the Galactic center region, von 
Linden et al. (1993) suggest that an accretion disk can 
also reproduce the basic structure seen in position- velocity 
plots. While this is certainly true for the inner region (x2- 
orbits are almost indistinguishable from circular orbits) it 
is not clear whether such a scenario can also explain the 
gas distribution and the different chemical properties in 
the entire CMZ. Based on the data presented in Papers 
I and II and this work, we will examine the questions 
related to large scale structure and dynamics closely in a 
forthcoming paper (von Linden et al. 1998). 

5.2. The origin of SiO 

Even though the SiO emission we observe is mostly as- 
sociated with cool gas, we consider it likely that it forms 
under conditions of high kinetic temperature and grain 
erosion by shocks (Martin-Pintado et al. 1992). Both the 
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Fig. 5. The chemical complexity of sources strong in SiO 
is demonstrated by our spectrum of M+0.83-0.18, cover- 
ing - 1.8 GHz and including both 29 SiO and 30 SiO. The 
spectrum shown is a composite of two single SEST low 
frequency resolution spectra. 

distribution discussed in the previous section and the ex- 
istence of the hot cloud core M+l.31-0.13 can be taken as 
evidence for this. From the fact that only one cloud in our 
sample, namely M+l.31-0.13, shows hot gas emitting in 
SiO lines, we can conclude that shocked gas forms dense, 
cooler cores fairly fast. Extrapolating from the analysis in 
Hollenbach (1988), we find a very short cooling time scale 
(•C 10 6 yr) for gas at a density of 10 4 cm~ 3 . A competing 
process is the condensation of gas phase SiO onto grain 
mantles. We estimate this 'freeze out' time scale to be of 
order 10 6 yr or slightly less (Rohlfs & Wilson 1996). Thus, 
gas with a high gas phase relative abundance of SiO has 
been cool for most of its lifetime. 

After more than 10 6 years, the gas slowly loses the 
chemical memory of having been shocked. From the typi- 
cal radial velocities of the clouds, we estimate that, during 
this time, these can move ~ 50 pc. If the very SiO-rich 
clouds are preferentially located on xi-orbits along the 
bar, these will remain sufficiently close in longitude to the 
region where they encountered the shock to be recognised 
as a distinct population. 

Since we find some SiO in all cloud cores, additional 
formation of gas phase SiO is required, probably by lo- 
cal turbulence and/or cloud-cloud collisions causing some 
shocks. This is a general characteristic of the cores of 
Galactic center GMCs. Ongoing local shock activity is 
likely to be also necessary as the main heating mechanism 
for the hot gas component seen in all clouds (see Flower 
et al. 1995 and the discussion in Hiittemeister et al. 1993). 

It is known that gas in post-shock regions, away from 
chemical equilibrium, is characterized by abundant com- 
plex molecules (e.g. Brown et al. 1988). This is exactly 
what can be observed toward the positions within the 
CMZ that are strong in SiO. We show just one example 
in Fig.||. In the Galactic disk, such complex spectra are 
typical for very confined regions, while in the CMZ they 
seem to be ubiquitous. This chemical complexity can only 
be maintained for ~ 10 6 years and thus requires frequent 
shocks. These may also be typical for the starburst envi- 
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ronment (see Henkel et al. 1987, Mauersberger & Henkel 
(1991), Mauersberger et al. 1991, Hiittemeister et al. 1997 
for the case of NGC 253). On larger timescales, an equilib- 
rium state is reached. This consists of a gas phase compo- 
nent of mostly diatomic molecules (Herbst & Leung 1989). 

Most time-dependent chemical models have been cal- 
culated for either cold, quiescent clouds (e.g. Bergin et al. 
1995) or 'classical' collapsing, star-forming hot cores (e.g. 
Brown et al. 1988, Caselli et al. 1993). Thus, detailed, 
model-supported recommendations of the molecules that 
should be used in future work to trace the chemical 
evolution of the the Galactic center cores cannot be 
given. CH3OH and SO2 appaer, however, to be promis- 
ing species. The former is a high temperature, high den- 
sity tracer requiring Tki n > 70 K to evaporate from dust 
grains. Its wealth of emission lines allows density and 
temperature determinations from multi-level studies. Ob- 
servations in external galaxies show that its abundance 
can differ widely even in the central regions of galaxies 
(Hiittemeister et al. 1997). SO2 is associated with shock 
chemistry and grain destruction. The abundance of the 
two molecules need not be correlated, as is shown by an 
evolutionary study of cores in the W3 region (Hclmich 
et al. 1994). Both molecules have lines that are strong 
enough to allow large scale mapping in the Galactic cen- 
ter region, and differences in their distribution would give 
further insight into the predominant processes operating 
on the Galactic center cores. 

6. Conclusions 

Based on our measurements and analysis of 33 molecular 
density peaks in the Galactic center in the J = 1 — > 0, 2 — > 
1 and 3 -> 2 transitions of C 18 0, the J = 2 -> 1 and 
5^4 transitions of 28 SiO and the J = 2 — > 1 transitions 
of 29 SiO and 30 SiO, we find: 

1. All sources are easily detected in all transitions of C 18 
searched for and in 28 SiO(2 — > 1), demonstrating that 
the properties of molecular peaks in the Galactic center 
region are markedly different from the Galactic disk, 
where thermal SiO emission is confined to very small 
regions in the vicinity of outflows associated with star 
formation. Local gas on the line-of-sight, distinguished 
from Galactic center gas by its narrow lines, is not 
detected in SiO. 

2. The rare isotopomcrs 29 SiO and 30 SiO are detected in 
all 12 studied sources. The J = 5 — > 4 transition of 
the main isotopomers is seen in 7 out of 8 sources. The 
line intensity ratio of 29 SiO and 30 SiO shows that the 
terrestrial isotope ratio, 29 Si/ 30 Si = 1.5, holds for the 
Galactic center region. 

3. From LVG model calculations applied to the C 18 
line intensity ratios ^2/1° an< ^ ^3/2°' co °l Pkin~ 
20 — 30 K) gas toward the molecular peaks has moder- 
ately high densities (n (H 2 ) ~ (1 — 4) 10 4 cm~ 3 ), while 



high kinetic temperatures of > 100 K correspond to H2 
densities that are an order of magnitude lower. This is 
contrary to what is found in the disk, where the cores 
of GMCs are usually hot, and is an indication that the 
'cool cores' in 'typical' Galactic center GMCs, away 
from Sgr A and SgrB2, are not presently forming high 
mass stars. 

4. Combining the results of LVG models for C 18 and 
SiO, and using C 18 as a tracer of total H 2 column 
density, a beam averaged SiO abundance is derived 
for all clouds. This varies significantly from source to 
source, ranging from ~ 0.3 10~ 9 to > 5 10~ 9 . The 
28 SiO(2 — > 1) transitions are optically thick, with to 
ranging from 2 — 8. 

5. Including information on the temperature structure of 
the clouds from NH3, it is shown that for most clouds 
a self-consistent solution accounting for the properties 
of C 18 and SiO is only possible if the bulk of the SiO 
emission arises in the cool, dense gas component. 

6. One source, M+1.31— 0.13, is different from all the oth- 
ers: SiO arises in the hot, thin gas component. Since 
this is also the cloud with the highest SiO abundance 
in the entire sample, the SiO formation process, proba- 
bly grain erosion by shocks, is likely to be still ongoing 
in this source. 

7. The SiO abundance in individual clouds is related to 
the large scale gas dynamics in the Galactic center re- 
gion. The highest abundances are found at Galactic 
longitudes of I > 0.8°, which can be identified with 
'collision regions' likely to encounter shocks in terms 
of the bar model of gas dynamics. 

8. As in the disk, SiO in the Galactic center region is 
likely to originate in shocks. Since the hot, thin post- 
shock gas forms dense, cool cores faster than SiO re- 
condenscs to dust grains, the SiO rich gas is cool for 
most of its lifetime. On a timescale of ~ 10 6 yr, the SiO 
molecules freeze out on grain mantles and the clouds 
lose their chemical memory of the shock. 
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